function printout_dispersion_eq_from_omega(m, CONSTS, plot_data)

    omega_p = CONSTS.omega_p;
    omega_H = CONSTS.omega_H;
    num_points = 1000;
    
    p_high = 10;
    num_intervals = 2000;
    use_combination = false;

    f_vec = linspace(omega_H/(2*pi), omega_p/(2*pi), num_points)';
    %determ_vec = zeros(size(f_vec));
    
    for i = 1:size(f_vec, 1)
        f = f_vec(i);
        LOCAL_CONSTS = CONSTS;
        LOCAL_CONSTS.f = f;
        LOCAL_CONSTS = update_CONSTS(LOCAL_CONSTS);

        z_det = fzeroall(@(p)imag(det(func_delta(p, m, LOCAL_CONSTS))), ...
        [1, p_high], num_intervals);

        %determ_vec (i) = det(y1); 
        
    end

    if(plot_data)
       
        figure; plot((2*pi)*f_vec, z_det); title('dispersion equation from \omega')
        
    end

end